Langevin approach to synchronization of hyperchaotic time-delay dynamics 
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In this paper, we characterize the synchronization phenomenon of hyperchaotic scalar non-linear 
delay dynamics in a fully-developed chaos regime. Our results rely on the observation that, in that 
regime, the stationary statistical properties of a class of hyperchaotic attractors can be reproduced 
with a linear Langevin equation, defined by replacing the non-linear delay force by a delta-correlated 
noise. Therefore, the synchronization phenomenon can be analytically characterized by a set of cou- 
pled Langevin equations. We apply this formalism to study anticipated synchronization dynamics 
subject to external noise fluctuations as well as for characterizing the effects of parameter mismatch 
in a hyperchaotic communication scheme. The same procedure is applied to second order differential 
delay equations associated to synchronization in electro-optical devices. In all cases, the departure 
with respect to perfect synchronization is measured through a similarity function. Numerical simu- 
lations in discrete maps associated to the hyperchaotic dynamics support the formalism. 

PACS numbers: 05.45.Xt, 05.45.Jn, 05.40.Ca, 05.45.Vx 
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I. INTRODUCTION 

In the last decades, synchronization of chaotic dynam- 
ics became a subject that has attracted of lot of attention 
0, HI, HI, IE IE IE 01 ■ I n fact , from a theoretical point of 
view, this phenomenon seems to contradict the inheriting 
sensibility to initial conditions of chaotic dynamics. On 
the other hand, the interest in this kind of phenomenon 
comes from the possibility of using the unpredictable 
chaotic trajectories as a carrier signal in communication 
channels. In this context, high dimensional systems with 
multiple positive Lyapunov exponents, i.e., hyperchaotic 
dynamics [i], have been proposed as a resource for im- 
proving the security in the communication schemes. Syn- 
chronization of hyperchaotic systems has therefore also 
become an area of active research [E [IE HE [IE Hi] . 

Chaotic dynamics described by differential delay equa- 
tions arise in the description of many different kind of 
situations, such as physi ology 114 
omy [T3| , laser physics |l8t Il9l. l2dl. 



bi olog y [161 ]. econ- 
2l|, mf, etc. As is 
well known, a high dimensional chaotic attractor char- 
acterizes these infinite dimensional systems. It has been 
shown that the Lyapunov dimension of the attractor is 
proportional to the characteristic delay time of the dy- 
namics HE HE HE Ull • Therefore, " synchronization of 
hyperchaotic nonlinear delay dynamics" has also been 
extensively explored from both a theoretical point of 
view as well as a resource for communication schemes 
HE HE HE HI]- A new aspect introduced in this case 
is the possibility of synchronizing two chaotic dynamics 
with a time shift, giving rise to the phe nomenon of an- 
ticipated HE HE Hi Ha HE Hi, HEM Hi, Hi, m (or 

retarded) synchronization, i.e., one of the chaotic systems 
(slave or receiver system) follows the chaotic trajectory 
of the other one (master or transmitter system) with an 
advanced (or retarded) time shift. 



In any real experimental setup where chaotic synchro- 
nization is observed one is naturally confronted with two 
undesirable effects that avoid reaching a perfectly syn- 
chronized regime. The characteristic parameters of both 
systems are not exactly the same, small parameter mis- 
match may induce clearly observable effects [IE Hi, Hi] . 
Also, departure with respect to the perfectly synchro- 
nized manifold may also be due to intrinsic noise fluc- 
tuations present in both systems [IE l45j - Both effects 
have been analyzed in the literature. Nevertheless, due 
to the chaotic character of the dynamics, in general it is 
hard to obtain an analytical estimation of these effects, 
which in fact may depend on the specific nature of the 
chaotic systems as well as on the coupling scheme used 
to achieve synchronization. 

The aim of this paper is to provide a simple analyti- 
cal description of the phenomenon of synchronization of 
hyperchaotic delay dynamics, considering realistic situa- 
tions such as the presence of parameter mismatch as well 
as the presence of intrinsic noise fluctuations in both syn- 
chronized systems. We demonstrate that this goal can be 
achieved when the synchronized dynamics are in a fully- 
developed chaos regime [46|. In this situation, the corre- 
sponding chaotic attractor does not have any stable peri- 
odic orbit and its basin of attraction fills out almost the 
whole available domain. These properties suggest that 
the nonlinear delay contribution terms of the chaotic dy- 
namics may be statistically equivalent to an ergodic noise 
source. As demonstrated in Ref. HE Hi] , this property, 
in a long time limit and depending on the characteris- 
tic parameter values, is in fact valid for a broad class 
of scalar delay dynamics. Therefore, our main idea con- 
sists in replacing the full set of coupled delay chaotic 
evolutions that lead to synchronization by a set of cor- 
related Langevin evolutions obtained from the original 
chaotic ones after replacing the nonlinear delay contribu- 
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tions by a noise term. Since the final Langevin equations 
are linear, their statistical properties can be obtained in 
an exact analytical way, providing a simple framework for 
characterizing the synchronization phenomenon. Depar- 
ture with respect to perfect synchronization is character- 
ized in terms of a similarity function |49j , which measures 
the degree of correlation between the quasi-synchronized 
dynamics. 

The paper is organized as follows. In Sec. II we review 
the conditions under which hyperchaotic delay dynamics 
can be represented by a Langevin dynamics. In Sec. Ill 
we analyze the phenomenon of anticipated synchroniza- 
tion when perturbed by external additive noises. In Sec. 
IV we analyze the effect of parameters mismatch in a 
hyperchaotic communication scheme [26]. In Sec. V we 
study a set of second order differential delay equations 
associated to an electro-optical laser device [43( under 
the effect of parameter mismatch and under the action 
of external additive noises. In all cases, we present nu- 
merical simulations that sustain our theoretical results. 
In Sec. VI we give the conclusions. 



II. LANGEVIN APPROACH TO 
HYPERCHAOTIC DELAY DYNAMICS IN THE 
FULLY-DEVELOPED CHAOS REGIME 

We will consider scalar nonlinear delay evolutions with 
the structure 

u(t) = -~,u(t)+(3f[u(t-T)]+l(t). (1) 

Here, 7 defines a dissipative time scale and T is the 
characteristic time delay. The parameter (3 controls the 
weights of the nonlinear function f(x), which is assumed 
to be oscillatory or at least exhibiting many different ex- 
trema 50]. The term I(t) represents an extra inhomoge- 
neous contribution that may corresponds to a stationary 
Gaussian white noise (Sec. Ill), a deterministic signal 
(Sec. IV) or even it may be a linear functional of the 
process u{t) (Sec. V). 

The parameter (3T may be considered as a complexity 
control parameter. In fact, when f3T 3> 1, the dynamic 
reaches a fully-developed chaos regime, where the sta- 
tionary statistical properties of the corresponding attrac- 
tor can be reproduced with a linear Langevin equation 
[47l |48| . This property can be understood by integrating 
Eq. {!]), in the long time limit (ft 3> 1), as 

u(t)» f dt'e-^ t ~ t ">f[f3u(t / -T) + I(f -T)], (2) 
Jo 

where u(i) = [u(t) — I(t)]//3, and l(t) is defined by 
I(i) = f Q t dt'e-' y( - t - t 'h(t'). Then, by writing the inte- 
gral operation over f(x) as a discrete sum, Jg dt'g(t') — > 
dt' g(jdt'), one realizes that u(t) may be considered 
as the result of the addition of many "statistically inde- 
pendent" contributions. This last property follows from 



the fact that f[f3x] oscillates so fast that it behaves as 
a driving random force. The characteristic correlation 
time [in units of time I/7] of the nonlinear force f[(3x] 
is of order 7//? [4?], |48|. When the correlation time is 
the small time-scale of the problem, i.e., much shorter 
than the characteristic delay time, 7//? <C 7T, and con- 
sistently much shorter than the characteristic dissipative 
time, 7//? <C 1, the statistical independence of the dif- 
ferent contributions follows. Notice that this last prop- 
erty is independent of the structure of the inhomogeneous 
term I(t), which only introduce a shift in the argument 
of /[/&]. 

From the previous analysis, in the parameters region 
[3T > 1 and j//3 <C 1, the central limit theorem [5l| tell 
us that u(t) is a Gaussian process. Due to this charac- 
teristic, this regime is also named Gaussian chaos. Only 
when I(t — T) is a non-linear function (functional) of u(t) 
[48| . the long time statistic may depart from a Gaussian 
one. 

The lack of correlation between the different contri- 
butions of the "chaotic force" f[f3u(t - T) + I(t - T)], 
allows to introduce the following Ansatz. In the fully- 
developed chaos regime, "the long time statistical prop- 
erties" of u(t) can be equivalently obtained from Eq. ([2]) 
after replacing the chaotic force by a noise contribution 
{Pf[u(t ~T) + I(t - T)] -► r)(t - T)), delivering the 
Langevin equation 

u(i) = -7ti(t)+ij(t-r) + I(i). (3) 

The noise rj(t) must have the same statistical properties 
as the delayed chaotic force. Since we are restricting our 
analysis to the Gaussian chaos regime-jt is sufficient to 
map the first two statistical moments [52J 

= lim /3/R0I, (4) 

t — *-oo 

V(t + r) V (t) = lim P 2 f{u(t + T)]f[u(t)}. (5) 

r — >oo 

The limit operation (limt_>oo) is introduced because the 
statistical mapping is only valid in the long time regime. 
The overbar symbol denotes an average over realizations 
obtained from Eq. {1} with different initial conditions. 
Equivalently, since the chaotic attractor is ergodic (by 
definition of fully-developed chaos regime) when the in- 
homogeneous term is statistically stationary [5lJ , this av- 
erage can also be considered as a time average. For ex- 
ample, limt—oo /[«(*)] = limf— oo(l/t) J dt'f[u(t% 

In agreement with the lack of correlation between the 
different contributions of the chaotic force, r)(t) can be 
approximated by a delta-correlated noise, 

^Mn-W)W)^AS(t-t'). (6) 

This white noise approximation applies for time inter- 
vals (t — f) larger than the characteristic time correlation 
(7//?) [in units of time I/7] of the chaotic force (47l. [48|. 

i.e.,7(t-*0>(7//J). 

The coefficient A measures the amplitude of f[u(t)]. 

Clearly, A must be proportional to (3 2 . Nevertheless, its 
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exact value is not universal and depends on the specific 
function f(x) [I?) as well as on the parameters that define 
the inhomogeneous term (see next Section). 

When 1(f) is defined by an external driving force, we 
will assume that 



[77(f) - v(t)W) = o, 



(7) 



i.e., the noise fluctuations representing the chaotic force 
and the inhomogeneous term are statistically uncorre- 
cted. The plausibility of this assumption follows from 
the fact that the white nature of r\ (f) only relies on the 
rapid oscillating nature of f[/3x) while it is not affected 
by the shift introduced by the functional 1(f). Under this 
condition, the inhomogeneous contribution only affects 
the mean value of the Gaussian profile associated to u(t). 

Finally, we will assume that 77(f) = 0. The validity of 
this condition only depends on the specific properties of 
the non-linear chaotic force. In fact, the rapid oscillating 
nature of f[f3x] allows to discarding the asymmetry intro- 
duced by 1(f) [Eq.(j2])]. As will become clear in the next 
section, the generalization to the case lim^oo (3f[u(t)] ^ 
can also be worked straightforwardly. 

With the noise correlation Eq. © , the stochastic evo- 
lution Eq. (J3j) becomes a (driven) Orstein-Ulenbeck pro- 
cess [51] . The statistical equivalence in the stationary 
and fully developed chaos regimes of the deterministic 
evolution Eq. (Q]) and the Langevin Eq. (without 
the inhomogeneous term) was proved in Refs. [47], H^j. 
Clearly, this stochastic representation does not provide 
any new information about the chaotic dynamics. Nev- 
ertheless, in the next sections we will use this equivalence 
for formulating a simple framework that allows us to get 
an analytical characterization of the chaotic synchroniza- 
tion phenomenon for different realistic circumstances. 



III. ANTICIPATED SYNCHRONIZATION 
SUBJECT TO EXTERNAL ADDITIVE NOISES 

In this section we will apply the previous Langevin 
representation of a hyperchaotic attractor to analyze the 
phenomenon of anticipated synchronization [30l l3ll |32| 
in the presence of external noise sources. We consider a 
complete replacement scheme (39| , defined by the coupled 
chaotic evolutions 

x(t) = - lx (t)+{3f[x(t-T)]+£ x (t), (8a) 

ij(t) = -7i/(<) +/?/[*(*)] +£„(*)• (8b) 

In this context, the variables x(t) and y(t) are referred 
as master and slave variables respectively. As before, 7 
is a constant dissipative rate, and (3 measures the weight 
of the delay-nonlinear contribution f(x). 

We have considered external noise sources, defined by 
the master and slave noises £ x (t) and £ y (t) respectively. 
We assume that their mean values are null (£ z (£)) = 0, 
where z = x,y, and (• • • ) denotes average over noise re- 
alizations. Furthermore, we assume that both noises are 



Gaussian, with correlations 

(t z (t)tAt'))=A zz ,5{t-t'). 



(9) 



The "diffusion" coefficients satisfy the positivity con- 
straint AxxAyy AxyAyx ^ 0, i,Axy — Ayx} |5ll] ■ These 

definitions allows to cover the case where both the mas- 
ter and slave variables are affected by intrinsic uncorre- 
cted fluctuations, A xy — 0, as well as the case of cor- 
related fluctuations A xy > 0. This last situation may 
be easily produced in any experimental setup. In par- 
ticular, for A xx — A yy — A X y, the noises that drive 
the master and slave dynamics are exactly the same, 
i.e., Eq. © with £ x (t) = £„(*). This property follows 
by diagonalizing the matrix of diffusion noise coefficients 

{{^ss , Axy } , {Axy , Ayy } } . 

As it is well known, in the absence of the external 
noises £ x (t) and £ y (t), the master-slave dynamics Eq. ([8]), 
after a time transient of order I/7, reach the synchro- 
nized manifold x(t + T) = y(t). Therefore, the slave vari- 
able anticipates the behavior of the master variable. The 
achievement of this state is clearly affected by the pres- 
ence of the external noises. The degree of departure with 
respect to the perfectly synchronized manifold can be 
measured with a similarity function, defined as [4^ | 



S(r) 



lim 

t — *oo 



{[x{t + r)-y{t)f) 

[<z 2 (u>(z/ 2 (u>] 1/2 



1/2 



(10) 



As before, the overbar denotes an average with respect to 
the system initialization or equivalently a time average 
in the asymptotic regime. In the absence of the external 
noises, this object satisfies S(T) = 0, indicating that the 
perfectly synchronized state x(t+T) = y(t) was achieved. 
In the presence of the noises, we expect S(T) > 0. 

The characterization of the behavior of the similarity 
function S(t) from the evolution Eq. ||5J) is in principle 
a highly non trivial task. The major complication come 
from the chaotic nature of the master and slave dynamics. 
Even in the absence of the external noises, in general 
it is impossible to get an analytical expression for the 
similarity function. Nevertheless, if both dynamics are 
in the fully-developed chaos regime, from the previous 
section, we know that a simpler representation may be 
achieved. The Langevin approach to Eq. ([5]) reads 

x(t) = -yx(t)+r](t-T)+Z x (t), (11a) 
y(t) = -7W(i) + + ( llb ) 

Notice that this equation corresponds to Eq. ([5]) with the 
replacement f[x] — > rj and maintaining the respective 
time arguments. As before, the effective noise rj(t) is 
defined by the correlation Eq. ([5]). Furthermore, we will 
assume that 77(f) — 0. We will deal the case, 77(f) 7^ at 
the end of this section. 

While the nature of Eqs. (fTTj) is completely different 
to that of Eqs. ©, in the asymptotic time regime these 
Langevin equations, without the external noises ^(f ) and 
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£ y (t), also reach a perfectly synchronized state. In fact, 
without the external noises, from Eq. (fTTj) it is possible 
to write (d/dt)[x(t) - y(t - T)] = -j[x(t) - y(t - T)], 
implying that after a time transient of order (1/7) the 
manifold x(t + T) = y(t) is reached. From the previous 
section we know that the statistical properties of the cor- 
responding master process Eq. (l8aj) and Eq. (jllap are the 
same. Then, we estimate the similarity function Eq. (|10|) 
associated to the nonlinear evolution Eq. ([8]) from the 
simpler linear Langevin evolutions Eq. (jlip . 

In the long time limit (-ft ^> 1), the master and 
slave Langevin evolutions can be integrated for each 
realization of the noise r](t), which represent the non- 
linear force, and external noises [£ x (i) and £ v (t)] as 

x{t) a / * dt'e-TC*"*') [r)(f ~ T) + £ x (*')], and as y(t) a 
J 'dt'e^( f -'')[j)(t') +€ y (t')], respectively. By using the 
correlations Eq. (J6]) [with rj(t) = 0] and Eq. (J9]), it follows 



lim (x 2 (t)) 



2 7 



lim (y 2 (t)) 

t — >oo 



27 



In a similar way, we get 
A 



lim 

t — >oo 



s(t + r)y(i)) = — exp[- 7 |r-T|]4 
27 

Therefore, the similarity function reads 
(1 + A **+*«v ) - (e-^- T l -f 



2 7 



(12) 



exp[-7|r|]. 



S(t) = V2 



(1 + ^*)V2(1 + 4fL)V2 

(13) 

This is one of the main results of this section. Notice that 
this expression only depends on one free parameter, i.e., 
the amplitude of the chaotic force A. On the other hand, 
this result relies in assuming the absence of any statis- 
tical correlation between the noise rj(t) representing the 
chaotic force and the external noises {£ x (t), £„(*)}• Con- 
sistently with the delta correlated nature of both con- 
tributions, Eqs. © and ©, we have also calculated the 
extra contributions to Eq. (|13[) that appear by assuming 
a delta cross correlation between both kind of objects. 
Nevertheless, the numerical simulations presented along 
the paper contradict the existence of any extra correla- 
tion, supporting the (previous) arguments that explain 
the statistical independence between the chaotic force 
and any external source, Eq. ([7]). 

From Eq. (|13p one can analyze different limits. In the 
absence of external noise we get 

vl/2 



1/2 



S{t) = V2(l - exp[- 7 |r - T\ 



(14) 



Note that this expression does not depends on the chaotic 
force amplitude A, it is defined only in terms of the 
local dissipation rate 7 and the delay T. Consistently, 
S(t) satisfies the anticipated synchronization condition 
S(T) = 0. When the external noise sources are taken in 
account, we get 

1/2 

V2 2 /"",_, A " e A '— . (15) 



S(T) 



(1 



This value measures the departure with respect to the 
perfectly synchronized manifold x(t + T) = y(t). Notice 
that the correlation A xy always decrease the value of 
S(T). This effect is exponentially diminished when in- 
creasing 7T. 



A. Anticipated synchronization of delay maps 

The previous results can be extended to the case of 
anticipated synchronization in discrete delay maps. We 
consider the maps that follow after discretizing the time 
variable in Eq. ([5]), t — nSt, and integrating both the 
master and slave evolutions up to first order in the dis- 
crete time step St. We get 



X n +1 



■bf(a 



y n+ x = ay n + bf(x n ) + 



(16a) 
(16b) 



Here, uq defines the characteristic delay step and the 
"dissipative" rate satisfies < a < 1. For each n, and 
are independent Gaussian distributed variables with 

(£n£m) = A zz 'Snm, (z = %,y), subject to the constraint 
A xx Ayy A X yAy X > 0, (A X y = Ay X ) [5 lj | . The parameters 
of the continuous time evolution Eq. ([8]) and the discrete 
map Eq. (fTBj) are related by 



P=j t , T = n St, A zz * = ^f. (17) 



l-o 

~sT' 



For the discrete map, the similarity function Eq. (|10[) is 
defined as 



S m = lim 

n — *-oo 



((Xn+m ~ Un) 2 ) 



1/2 



(18) 



As we will show in the next examples, this object can 
be fitted by Eq. JT3]) with the mapping Eq. (JTTJ) . This 
property is valid when the discrete map provides a good 
approximation to the continuous time evolution. In fact, 
when awl, i.e., jSt <C 1, the map Eq. (| 16[) can be read 
as a numerical algorithm for simulating the continuous 
time evolution Eq. (TlTj) . Nevertheless, we remark that 
Eq. (fT(3|) can also be analyzed without appealing to the 
parameter mapping Eq. (|17[) . i.e., the Langevin approach 
can also be formulated for discrete time dynamic. S m can 
be estimated after replacing the chaotic force in Eq. (TI| 
by a discrete noise, f(x) — > r) n , with rfc = and rf^jT^ 
AS nm . We get 



A + Ai 



n-^oo I — a 



lim {yl) = (19) 



The master-slave correlation reads 



lim ((x n+m y n )) = -{Aa} m ~ n °\ + A xy a^). (20) 

n— too 1 — a z 
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Then, the (discrete) similarity function reads 



(1 



2A 



vv ) ( fl |m— no| _|_ A xy ^Im!) 



(1 + ^)1/2(1 + ^)1/2 



1/2 



(21) 



Consistently, by using Eq. (TT7|) and the mapping between 
the chaotic force amplitudes 



(22) 



A = ASr 



from Eq. (|2Tj) to first order in one recovers Eq. (fl"3|) . 
Furthermore, the conditions that guarantee the validity 
of the Langevin representation in the continuous time 
case, i.e., f3T 1 and 7//? <C 1, from the mapping Eqs. 
(fTF]) . here read bn > 1 and (1 - a)/b < 1. 



B. Numerical results 

To check the validity of the previous results, we con- 
sider an Ikeda-like delay-differential equation [lHOjj], i.e. 
Eq. (H|) with f[x] = s'm(x). Its associated discrete map 
[Eq. (HI] read 

x n+ i = ax n + bsin(x n - no ) + (23a) 
Vn+i = ay n + 6sin(x n ) + (23b) 

To obtain the following results, we generate a set of re- 
alizations from the map Eq. (|23|) by considering differ- 
ent random initials conditions in the interval (— 7r,7r) for 
both, the master and the slave variables. By averaging 
over this set of realizations (« 10 4 ), we determine nu- 
merically the similarity function Eq. (fTS)) . To check the 
ergodic property of the corresponding chaotic attractor, 
we repeated the numerical calculations by averaging over 
time a single trajectory with an arbitrary set of initial 
conditions. Consistently, we obtained the same results 
and characteristic behaviors. 

In Fig. 1 we show the similarity function when the ex- 
ternal noises are absent, i.e., the similarity corresponding 
to the deterministic map. In agreement with this condi- 
tion, we notice that S no = 0, implying that the manifold 
Xm+no — Dm characterizes the asymptotic behavior. Fur- 
thermore, we find that both the expression for the con- 
tinuous time case [Eq. (fT3")) with the mapping Eq. fTF])]. 
as well as the similarity function of the map [Eq. (f2"Tj) ] 
are indistinguishable from each other (in the scale of the 
graphic) and both correctly fit the numerical behavior. 
In the inset, we show the stationary probability distri- 
bution of the master process x n . In agreement with our 
considerations, this distribution can be fit with a Gaus- 
sian distribution. From its width, and by using Eq. (|19J) 
[or Eq. (fT2"|) ] we estimated the value of the chaotic force 
amplitude, A ss 3.5. 

In the absence of the noises, we corroborate that by 
taking the function f[x] = sm(x + <f), independently of 
the value of the phase <p 1 the same statistical behaviors 
follow. This result confirms the arguments presented in 
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FIG. 1: Similarity function of the coupled delay maps Eq. (|23[) 
without noise, A xx — A yy — A xy — 0. The parameters are, 
no = 15, b — 3, and a = 0.975 (squares) and a = 0.9 (cir- 
cles). The lines correspond to the fitting Eq. (|14[) joint with 
the mapping Eq. (|17|l [with St — 1], indistinguishable from 
Eq. (|21|l . The inset corresponds to the stationary probability 
distribution P(x) of the process x n (with a — 0.9). 



05 




1 1 r- 

15 30 45 60 



FIG. 2: Similarity function of the coupled delay maps Eq. (123[) 
subject to external Gaussian noise fluctuations. The param- 
eters are no = 15, a — 0.9, and 6 = 3. The noise parameters 
are A xx ~ A yy — 0.25 and A xy — 0. The line corresponds 
to the fitting Eq. (|14[) joint with the mapping Eq. (|17[1 [with 
5t = 1], indistinguishable from Eq. (|2ip . The chaotic noise 
amplitude results A w 3.5. The inset corresponds to a master 
(x n , full line) and slave (y n , dotted line) realization. 



the previous section [Eq. ©] about the statistical invari- 
ance of the chaotic force under a shift of its argument. 

In Fig. 2 we show the similarity function when the 
master and slave dynamics are affected by two uncor- 
related external noises. As expected, we found that 
S no > (no = 15). In the inset, we show a character- 
istic master-slave realization. Even in the presence of 
the external noises, both trajectories are approximately 
the same (the slave anticipate the master trajectory). In 
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FIG. 3: Similarity function of the coupled delay maps Eq. (|23[) 
driven by different noises. The circles correspond to A xx = 
A yy — 3.5 and A xy — 0, while the squares correspond to 
A xx = Ayy — A xy = 3.5. The parameters of the map are 
the same than in Fig. 2. The lines correspond to the fitting 
Eq. (fl4|) joint with the mapping Eq. dT7)l [with 5t = 1], in- 
distinguishable from Eq. (|21|l . The chaotic noise amplitude 
results A m 4.3. 



contrast with the previous figure, here the fitting to the 
similarity function depends explicitly on the chaotic force 
amplitude A. We found that the value of A that provides 
the best fitting is consistent with the one found for the 
deterministic map (Fig. 1), i.e., A w 3.5. Then, in this 
case, the inequality A {A xx , A yy } is satisfied, implying 
that the fluctuations induced by the determinist chaotic 
dynamic are much larger than those induced by the ex- 
ternal noise sources. 

Maintaining all the parameter values corresponding to 
Fig. 2, we analyzed the case A xx = A yy = A xy . In 
this situation, the noises that drive the master and slave 
dynamics are exactly the same, i.e., Eq. (|23p with ^ = 
£^ (as in the continuous time case, this property follows 
by diagonalizing the matrix of diffusion noise coefficients 
{{A xx ,A xy }, {A xy , A yy }}). We found that the similarity 
function is almost indistinguishable from that of Fig. 2. 
Independently of the external noise correlations, in both 
cases our approach provides a very good fitting of the 
numerical results. 

In Fig. 3, by maintaining the parameters of the de- 
terministic map, i.e., (a, b, no), we increased the ampli- 
tude of the external noises, such that {A xx , A yy } w A. 
Then, the external noise-induced fluctuations are of the 
same order as the intrinsic chaotic dynamical fluctua- 
tions. Even in this limit, our approach provides an ex- 
cellent fitting of the numerical results. Both, the case of 
uncorrelated noises (A xy = with A xx — A yy ) and the 
case of completely correlated noises (A xx = A yy = A xy ) 
are considered. In both situations, since the external 
noise amplitudes are larger than in Fig. 2, the value of 
S no increases, indicating a weaker (anticipated) synchro- 
nization between the master and slave dynamics. 



In the case of completely correlated noises, A xx = 
Ayy = A^y, the similarity function develop two (local) 
minima, one at m — no and other at to = 0. This feature 
follows from the competition between two different syn- 
chronizing mechanisms, induced by the chaotic dynamics 
and the external noises, respectively. In fact, as the noise 
that drives the master and slave dynamics is the same, 
its action tends to synchronize both dynamics without 
any time shift [see Eq. (j23f with £jj = producing the 
dip at to = 0. In agreement with this argument, from 
Eq. (f2"Tj) [or equivalently Eq. (|13[) ] one can deduce that 
only one minimum at m = will be appreciable in the 
similarity function when A <C A xx — A yy = A xy , i.e., in 
the limit of high noise intensity. 

As in Fig. 2, the similarity function depends on the 
amplitude A of the chaotic force. Here, the best fitting 
to the similarity function is obtained with A ss 4.3. This 
value is larger than the one associated to the determinis- 
tic dynamics [Fig. 1] or the case of weak external noises 
[Fig. 2]. Then, in general it is necessary to consider that 
the chaotic force amplitude A is also a function of the ex- 
ternal noise intensities, A = A[{f(x),b},A xx ,A yy )]. We 
remark that in concordance with the Langevin represen- 
tation, the value of A determined from Eqs. (flU]) to (|2"Tj) 
is the same. By using the parameters of Fig. 3, we found 
a moderate dependence on the external noise intensities, 
i.e., the maximal variation of A with the amplitude of 
the noises does not exceeds thirty percent (30%) of the 
deterministic map value (Fig. 1, A « 3.5). The depen- 
dence is smooth but non-monotonous. We found that A 
saturates to a fixed value (A » 4.5) when increasing the 
external noise amplitudes, A < {A xx , A yy }. 



C. Chaotic force with a non-null average value 

Our previous theoretical calculations rely on the as- 
sumption that the average (over realizations or its sta- 
tionary time average) of the chaotic force is zero, i.e., 
f[x(t)] = 0. When the dynamic, that does not include 
the non-linear contribution, is purely dissipative [Eq. ([1]) 
with f(x) — > 0], the validity of this assumption requires 
that f(x) takes symmetrically positive and negative val- 
ues. Functions that do not satis fy th is property also arise 
in real experimental setups [2ll. ]2rj |43|. In these cases, 
due to the pure dissipative nature of the evolution Eq. ([8]) 
[or Eq. (fl"6|) ]. in the long time limit the master-slave re- 
alizations will fluctuate around a non-null fixed value. 

This situation can be managed by writing f(x) = 
[,f(x) — f{x)} + f(x). Then, the previous theoretical cal- 
culations can be easily extended by replacing [f(x) — 
f(x)] — > r](t) [with n(t) = 0] and by maintaining the 
extra contribution f(x) in the final Langevin representa- 
tion. For example, taking the functions f(x) = sin 2 (:r) or 
f(x) = cos 2 (a;) in Eq. ©, from the Langevin represen- 
tation Eq. (fTTj) . it is possible to deduce that in the fully- 
developed chaos regime the master-slave realizations will 
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fluctuate around f3/(2j). 

When the dynamic without the non-linear contribution 
can by itself induce symmetric dynamical oscillations, the 
symmetry requirement on the function f(x) may be elim- 
inated [see Sec. V]. 



IV. EFFECT OF PARAMETERS MISMATCH IN 
A HYPERCHAOTIC COMMUNICATION 
SCHEME 

Hyperchaotic delay dynamics may be used as a re- 
source for secure encoded communication. Different 
schemes have been proposed, in all the cases, it is argued 
that the security of the communication channel may be 
improved by increasing the dimension of the hyperchaotic 
attractor. Here, we study the dynamics [26j 

x{t) = -jx(t)+{3f[x(t-T)]+Wl(t), (24a) 
y(t) = -iy{t)+P'f[x{t-T')]. (24b) 

The variable x(t) is the carrier signal where the mes- 
sage is encoded. The external feed Wl(t) is defined by 
DJl(t) = [dm{t)/dt + jm(t)], where m(t) is the message 
to be transmitted. The variable y(t) is the receiver. It 
is easy to demonstrate that in the stationary regime the 
state x(t) — y(t) = m(t) is reached, implying that the 
receiver is able to unmask the message encoded in the 
hyperchaotic dynamics of x(t). This condition is only 
satisfied when the characteristic parameters of the trans- 
mitter and receiver evolutions [Eq. (12"4"| ] are the same, 
i.e., 7' = 7, = (3, and T' — T |53|. In any real ex- 
perimental situation, it is expected that this condition 
is not fulfilled, i.e., the decodification of the message is 
performed in the presence of an unavoidable (finite) pa- 
rameter mismatching. For simplicity, in the present sec- 
tion we will not consider the action of any external noise 
perturbation source. 

In order to achieve a general characterization of the 
influence of the parameters mismatch, we take Wl(t) = 0. 
As in the previous section, we will use the similarity func- 
tion S(t) [Eq. (fTU)l] as a measure of the degree of syn- 
chronization between the emitter and receiver variables. 

When the transmitter dynamics is in the fully- 
developed chaos regime, we can extend the Langevin ap- 
proach to the present situation. In order to estimate the 
similarity function, we replace the chaotic coupled evo- 
lution Eq. (|2"4|) by the Langevin equations 



x(t) 

m 



-~,x(t)+r)(t-T), 
-</»(*) + V^v(t - T'), 



where the parameter a is defined by 

t 

ft 



(25a) 
(25b) 



(26) 



As before, the noise rj(t) is characterized by rj(t) 
V(t)v(f) = A6(t - t') [Eq. % 



0, and 

The similarity function 



S{t) associated to Eq. (|25[) . can be determine analytically 
by a straightforward calculation. We get 



S(r) 




y ( 7 + y) 



-exp[-7|r 



A| 



1/2 



(27) 



where we have defined the delay time mismatch 

A = T - T', (28) 
and the rate 7 is defined as 
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7 for t > A 
7' for t < A 



(29) 



We notice that due to the normalization constants in the 
definition of the similarity function Eq. (fT0|) . the final 
expression Eq. (|27|) does not depend explicitly on the 
chaotic force amplitude A. The similarity function only 
satisfies 5(0) = 0, when the transmitter-receiver parame- 
ters are exactly the same. In this situation, the manifold 
x{t) = y(t) characterize the stationary regime. 



A. Parameter mismatch in discrete delay maps 

The previous results can be extended to the coupled 
maps obtained from the evolution Eq. (|24p after a first 
order Euler integration, i.e., 

x n+ i = ax n +bf(x n - no ), (30a) 
y n+ i = a'y n + b'f(x n _ n > o ). (30b) 

Here, the transmitter and receiver parameters of both, 
the continuous and the discrete time evolutions, must be 
related as 

1-a b 
7 = — — , P 
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St 1 
1-a' 

St 



St' 



T = n$5t, 



(31a) 







5t' 



T' = n' St, (31b) 



where St is the characteristic discretization time step. For 
the discrete maps, the similarity function Eq. (|18p reads 



/ (1-a' 2 ) 
a(l-a 2 ) 



+ 



/a(l-a 2 ) 



(1 - a r 



2V(l-a 2 )(l-q' 2 ) , mn | 
(1 - aa') 



1/2 



where we have defined 



m = n - n , 



and the "dissipative rate" 



b' 



_ j a for m > mo 
a' for to < too 



(32) 



(33) 



(34) 



Notice that under the mapping Eq. (|3"Tj) , the definition of 
the parameter a in Eq. (|33p is consistent with Eq. 
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[ft 
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m 



FIG. 4: Similarity function of the coupled delay maps Eq. (|35[) 
(circles). The parameters are no = 15, a = 0.9, 6 = 3 and 
n' = 20, a' — 0.8, b' = 5. The line corresponds to the fit- 
ting Eq. p7|l joint with the mapping Eq. pip [with St — 1], 
indistinguishable from Eq. (I32|l . The inset corresponds to a 
realization of x n (full line), and y n (dotted line). 



B. Numerical results 



Here we consider the coupled chaotic maps 



X rL J r \ — CLX n -\- h COs(Xro— no ) : 

y n+ i = a'y n + b' cos(x n - n > a ). 



(35a) 
(35b) 



In the inset of Fig. 4, we show a characteristic realiza- 
tion of the emitter and receiver variables. By averaging 
over different realizations (~ 10 4 ), we get the similar- 
ity function. We notice that due to the specific values 
of the characteristic parameters, the receiver "synchro- 
nizes" with the "past" of the transmitter variable. In 
fact, S m attains it minimal value at tuq = — 5. Further- 
more, due to the difference in the parameters b and &', 
the amplitude of the receiver fluctuations are bigger than 
those of the transmitter. 

In contrast with Fig. 2, here the similarity function is 
not symmetrical around its minimal value. This asym- 
metry has its origin in the mismatch between the "dissi- 
pative rate" parameters a and a' . 

Clearly, the Langevin approach provides a very well 
fitting to the numerical simulations. Furthermore, it al- 
lows to characterize the influence of the parameters mis- 
match on the synchronization of the transmitter and re- 
ceiver variables. The analytical quantification of this ef- 
fect can be obtained from the similarity function Eq. (|2"T|) 
[or Eq. ([32)) ] evaluated at the origin 



S(0) 




,1/2 



(7 + 7') 



exp[-7|A| 



By a direct inspection of this expression, we realize that 
the dependence of S(0) on the prime parameters is non- 



monotonous, S(0) develops different minima when vary- 
ing the receiver parameters. In agreement with the re- 
sults of Ref. [43|, assuming that it is possible to adjust a 
given receiver parameter, the previous expression allows 
us to choose the best value that maximizes the synchro- 
nization phenomenon. 



V. SYNCHRONIZATION OF SECOND-ORDER 
NON-LINEAR DELAY DIFFERENTIAL 
EQUATIONS 

In the previous sections, we analyzed the phenomenon 
of chaotic synchronization for dynamics generated by first 
order delay equations. Nevertheless, second order delay 
differential equations also arise in the in the description of 
real experimental setups. By second order we mean equa- 
tions whose linear dynamical contributions are equivalent 
to second order (time) derivative evolutions. Here, we 
demonstrate that the Langevin approach also works in 
that case. 

Following Ref. [43[ we consider the evolution 

z(f) + — + ^ [ x(s)ds = 0cos 2 [x{t-T)]+i x (t), (37) 

y(t)+^T-+w (\{.B)ds = (3' C os 2 [ X (t-T')]+Zy(tW8) 
7 Jo 

which describe synchronization in a set of coupled electro- 
optical laser devices. The integral contributions propor- 
tional to (1/0) and (1/9') indicate that the linear evolu- 
tions defined by the left hand side of Eqs. ([57)1 and ([55)1 
are equivalent to a set of second order derivative differen- 
tial equations. In addition to the parameters mismatch, 
we also consider the action of external additive noises 
£ x (f) and £, y (t), whose mutual and self-correlations are 
defined by Eq. In order to simplify the final expres- 
sion here we not consider any mismatch in the phase of 
the chaotic forces (53j . 

Under the same conditions than in the previous sec- 
tions, we assume that x(t) and y(t), in the long time 
limit, attain the fully-developed chaos regime, allowing 
us to replace the nonlinear delay forces by noise contri- 
butions with the same time arguments, delivering 

x(t) 1 /"* 

+ — + - / *( s ) ds = v(t-T)+£ x (t), (39a) 
7 9 J Q 

y(t) + ^T- + W fv{s)ds = V^v(t-T')H y (t), (39b) 
7 6 Jo 

where as before, a = (j3 /(3) 2 - The correlation of r](t) is 
again defined by Eq. ([6]), and we take 77(f) = 0. Notice 
that in spite that the nonlinear contribution cos 2 [x] is al- 
ways positive, when {9, 9'} < 00, the second order linear 
evolution introduces self dynamical oscillations that im- 
ply that its effective action averaged over realizations (or 
its stationary time average) must be taken as zero. This 
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property breaks down when the integral contributions are The similarity function [Eq. (|10[) ] reads 
absent, i.e., in the limit 6 = 6' = oo. 

The Green functions associated to the linear evolutions 
Eq. (f3"9"| can be obtained straightforwardly in the Laplace 
domain, being defined by the addition of two exponential 
functions. After integrating formally both equations for 
each realization of the noises, it follow 



lim(x 2 (i)> = i(A + A xx ), (40a) 
lim(p(i)) = ^(aA + A yy ). (40b) 



S(r) 









-4/i 


h(t - A) ^ 








) 1/2 (i 







1/2 



(41) 



Here, A = T — T' . The function S(t) is defined as 



2,0 



S(t) 



= p -57l*l 



cosh($|t|/2)-i>2sinh($|i|/2) , (42) 



4> 



where the dissipative rate 7 is defined by 



7 i/ t > 
7' i/ i < 



For t > 0, the frequency $ reads 



* = 7\/l 



6»7 : 



*/ t>0, 



(43) 



(44) 



while the dimensionless parameter ^ is 
2 



P = 1 



+ 6') V7 77 ' 



i/ t > 0. (45) 



For < < 0, both $ and f are defined by interchanging 
6^6' and 7 <-» 7' in the previous expressions. Finally, 
in Eq. (|4"Tj) . we have also defined the dimensionless pa- 
rameter 



r^i(6- 



{6H6'i) + [{6 + 0')(7 + 7') - 2] + (0'7/y0) 



which is symmetric in the emitter and receiver parame- 
ters. 

In order to check the validity of the previous approach, 
we consider the discrete maps associated to Eqs. (f3"T|) and 
by Euler integration 

n 

x n+ i = ax n - ^22 x i + bcos2 ( x n-n ) + (47a) 

1=0 

n 

y n+1 = a'y n -cj'22 Vj + b' cos 2 + (47b) 

1=0 



0,0 




-5 5 



1 — i — < — i — < — i — • r 

100 -75 -50 -25 



m 



I 1 I 1 I f 

25 50 75 100 



FIG. 5: Similarity function of the coupled delay maps Eq. (|47[) 
(circles) without noises, A xx = A yy = = and without 
mismatching, n' = no, a' = a, b' = b and a/ = w. The 
parameters are no = 20, a = 0.9, 6 = 3 and u) = 0.05. The 
line corresponds to the fitting Eq. (|41|l under the mapping 
Eq. (|48p [with St = 1] . The inset corresponds to the stationary 
probability distribution of the process x n . 



(46) The (transmitter) parameters mapping reads 



7 



St 



6 = St 



(1-a) 







(1-a)' 



(48) 



The same relations are valid for the prime (receiver) pa- 
rameters. The noises and joint with the corre- 
sponding mapping for their amplitudes are defined below 
Eq. (fTB")) . The chaotic force amplitude of the continuous 
and discrete time evolutions are related by A = A/ St. 

In Fig. 5 we plot the similarity function [Eq. (fT8")) ] as- 
sociated to the coupled maps Eq. (|4"T|) in absence of pa- 
rameters mismatch and without the external noises. The 
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T 1 1 1 . 1 . 1 . 1 . 1 1 1 r 
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0,4- m 



1 — | — i 1 — i — i — i — i — i — | — i — | — i — | — i — 
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m 

FIG. 6: Similarity function of the coupled delay maps Eq. (147[) 
(circles) with parameters mismatch and under the action of 
external noises. The parameters are no = 20, a — 0.98, 6 = 3, 
u = 0.01, and n' = 22, a' = 0.95, b' = 4, J = 0.0125. The 
noise parameters are A xx = A yy — 0.0225 and A xy — 0. The 
line corresponds to the fitting Eq. (|41[) joint with the mapping 
Eq. (1481) [with St — 1] . The inset corresponds to the similarity 
function in absence of noises and without mismatching. 



analytical result Eq. (|4Tj) . under the mapping Eq. ([35]). 
provides an excellent fitting of the numerical results. In 
contrast with the previous sections (see for example Fig. 
1), here the similarity function develops an oscillatory be- 
havior, its origin can be associated to the integral contri- 
butions in the dissipative dynamics of Eqs. (f37|) and (|38|) . 
Their characteristic frequency follows from Eq. (l4"4"|) . 

In order to check the achievement of the fully- 
developed chaos regime [s3 |. in the inset we show the 
stationary probability distribution of x n . Consistently, 
this distribution can be fitted with a Gaussian distribu- 
tion. 

In Fig. 6, for other set of characteristic parameters, 
we plot the similarity function in presence of parameter 
mismatch and external noise sources. As in the previous 
case, the fitting Eq. (|41 [) correctly match the numerical 
results. The effective chaotic force amplitude is A w 1.2. 
The asymmetry of S m around its minimum value follows 
from the parameter mismatch between the evolution of 
x n and y n . 

In the inset, maintaining the parameters correspond- 
ing to the evolution of x n , we plot the similarity function 
in absence of both, parameters mismatch and the exter- 
nal noise contributions. In this case, the chaotic force 
amplitude [determine from Eq. (|4TJ| ] is Am 1. 

As in the previous section, a remarkable aspect of our 
analytical results is that they allows us to know the in- 
fluence of different realistic effects on the chaotic syn- 
chronized manifold. For example, for a given set of fixed 
parameters, which may include the amplitude of the ex- 
ternal noises, one can determine the value of the rest of 
the parameters that minimizes the similarity function at 
t = 0, giving rise to a maximization in the degree of syn- 



chronization between the emitter and receiver variables. 
This kind of analysis follows straightforwardly from our 
analytical results. 



VI. CONCLUSIONS 

In this paper, we have characterized the phenomenon 
of chaotic synchronization in scalar coupled non-linear 
time delay dynamics. Our formalism relies in recogniz- 
ing that in a fully-developed chaos regime, the trajecto- 
ries associated to a broad class of hyperchaotic attrac- 
tors are statistically equivalent to the realizations of a 
linear Langevin equation. This equivalence can be estab- 
lished when the function that defines the driven chaotic 
delay force is an oscillatory one, such that its dynami- 
cal action can be represented by a delta-correlated noise. 
Given this Langevin representation, the coupled nonlin- 
ear delay evolutions, where the phenomenon of chaotic 
synchronization happens, are replaced by a set of lin- 
ear stochastic equations where the noise that represents 
the chaotic force maintains the corresponding time argu- 
ments. The statistical properties of the Langevin equa- 
tions can be obtained analytically, providing an excellent 
estimation of the stationary statistical properties of the 
synchronized manifold. 

Using the Langevin representation, we analyzed the 
phenomenon of anticipated synchronization in the pres- 
ence of external additive noises. We also characterized 
the effect of parameters mismatch in a hyperchaotic com- 
munication scheme. Second-order delay equations associ- 
ated to an electro optical device were also characterized. 
The analytical predictions of the Langevin approach cor- 
rectly fit numerical simulations in discrete coupled non- 
linear delay maps associated to the corresponding con- 
tinuous time evolutions. 

In all cases, the degree of synchronization (between the 
master-slave or emitter-receiver variables) was character- 
ized through a similarity function, defined in terms of the 
correlation between the synchronizing systems. When 
the departure from perfect synchronization is due to a 
parameter mismatching, the fitting to the similarity func- 
tion does not involve any free parameter. When the ac- 
tion of external noises is considered the fitting depends 
on the effective chaotic force amplitude. 

Our results are interesting from both, theoretical and 
experimental point of view. From our analytical expres- 
sions it is possible to evaluate under which conditions a 
given undesired effect can be minimized by controlling 
the rest of the parameters. On the other hand, in the 
context of hyperchaotic communication schemes, while 
high dimensional systems increase the complexity of the 
masking signals, our results show that the corresponding 
statistical properties may adopt a simple analytical form. 
In fact, by measuring the similarity function our results 
allow to infer the value of some of the characteristic pa- 
rameters of the hyperchaotic delay dynamics. 

The present study may be continued in different rel- 
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evant directions such as the extension of the Langevin 
representation beyond the fully-developed chaos regime 
(non-Gaussian chaos) as well as for non-scalar chaotic 
dynamics. Furthermore, the characterization of the de- 
pendence of the effective chaotic force amplitude with the 
external noise parameters is an open interesting issue. 
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